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THE  GROWTH  OF  HOT  REDUCED  DENSITY  CHANNELS 
IN  GASES  DUE  TO  TURBULENCE  AND  HEAT  CONDUCTION 


I.  INTRODUCTION 

A  hot  channel,  formed  by  either  laser  or  ohmic  heating  of  a  gas,1  2  expands  initially  due  to  over¬ 
pressure.  If  spatial  asymmetries  exist  in  the  initial  heating,  vorticity  is  created  and  degenerates  into  tur¬ 
bulence  throughout  the  channel3.  A  second  phase  of  growth  then  follows  the  attainment  of  pressure 
equilibrium  as  outside  gas  mixes  into  the  hot  channel.  Viscous  effects  cause  the  turbulence  to  decay 
until  heat  conduction  becomes  significant.  As  the  gas  responds  to  this  heat  flow  a  final  growth  and 
filling-in  of  the  channel  occurs. 

We  predict  the  scaling  of  the  channel  radius  with  time  during  both  later  phases  of  growth.  To  do 
this  we  use  the  governing  equations,  the  assumption  of  pseudo-self-similar  expansion  and  in  the  first 
case  a  simple  "mixing  length"  model  of  the  turbulence. 

II.  GROWTH  DUE  TO  TURBULENCE 

In  turbulent  flow  a  pattern  of  chaotic  eddies  with  a  spectrum  of  sizes  is  superimposed  on  the 
mean  fluid  motion.  If  the  fluid  velocity  is  represented  as  a  mean  value,  vr,  plus  a  relative  fluctuation, 
v'(,  such  that 

Vr  “  Vr  +  V'r. 

then  the  rms  value  of  the  fluctuation  velocity  is  a  measure  of  the  intensity  of  the  turbulence 

v,  -  (v>2. 

Turbulence  embedded  in  a  uniformly  moving  fluid  will  rapidly  become  isotropic 

vr  -  v,  -  v.  -  v. 
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We  shall  model  the  channel  as  entirely  filled  with  isotropic  turbulence  and  undergoing  a  cylindrical^ 
symmetric  average  expansion. 

Turbulence  transports  fluid  properties  by  small  scale  convection.  However,  this  effect  may  be 
modeled  as  a  large  scale  version  of  transport  at  the  molecular  level  and  treated  as  an  anomalous 
diffusion.  We  suppose  small  volumes  of  fluid  move  a  distance  /  (Prandtl  eddy  length)  before  breaking 
up  and  losing  their  identities.  The  product  of  this  "mean  free  path"  and  the  relative  velocity  is  known 
as  the  eddy  diffusivity  and  becomes  the  effective  mass  diffusivity,  kinematic  viscosity,  or  thermal 
diffusivity 

D  ”  v  —  a  “  /v. 

(Slightly  different  values  of  /  actually  apply  to  each  situation.) 

Experiments  show  isotropic  turbulence  has  particularly  simple  properties.'*  If  a  uniform  flow 
passes  through  a  grid,  the  resulting  turbulence  becomes  isotropic  a  short  distance  downstream  from  the 
grid.  The  intensity  of  the  turbulence  decays  beyond  this  virtual  origin  according  to 

1  1 

where  the  second  step  follows  from  the  constant  mean  velocity.  The  eddy  length  however  grows 

/  —  Vx  —  v7. 

We  may  therefore  ascribe  constant  transport  properties  to  the  medium  as  long  as  the  turbulence  per¬ 
sists.  The  value  of  this  constant  depends  however  on  the  initial  conditions  under  which  the  turbulence 
was  created,  in  this  case  on  the  grid  wire  size. 

We  may  take  spatial  averages  of  the  variables 

p-±fpdY  (I) 

v--pj fvrfP  (2) 

T  “  ~~  J*  p  TdV  (3) 

n  V  *v 
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where  we  select  the  averaging  volume.  V,  as  larger  than  the  scale  of  the  turbulence  but  smaller  than  the 
channel  scale.  The  form  of  the  governing  equations,5  when  expressed  in  terms  of  the  average  proper¬ 
ties,  is  found  by  application  of  Reynolds  procedure6.  For  example  we  substitute  in  the  continuity  equa¬ 
tion  as  follows 

T~(p  +  p  )  +  V-[(v  +  v’)  (p  +  p')l  *  0  (4) 

ot 

and  then  average  the  entire  equation  resulting  in 

+  V-(pv)  +  V-(£V)  -  0.  (5) 

ot 

In  recognition  of  the  transport  effect  of  the  turbulence  (and  following  Boussinesq’s  momentum 
transport  treatment7)  we  make  the  following  replacement  for  the  correlation  term 

pV  ”  -aVp  (6) 

with  the  result  that 


+  Vpv  —  aV^p  =  0. 
ot 

(7) 

The  conservation  of  energy  equation  may  be  treated  similarly  and  the  correlation  terms  replaced  by  a 

diffusive  one.  The  governing  equations  thus  become 

+  V-p  v  -aV2p  -  0  conservation  of  mass 

Ot 

(8) 

PCy 

~  +  v  •  V  fj  +  p0RT0V  v  -  V  -  |pcy«  V  7j  -  0  conservation  of  energy 

(9) 

P T  •  p„T0  equation  of  state 

(10) 

This  set  of  equations  incorporates  a  number  of  simplifications.  The  Navier-Stokes  equation  (conserva¬ 
tion  of  momentum)  has  been  eliminated  by  assuming  constant  pressure.  A  symmetric  expansion 
creates  no  new  turbulence  and  the  energy  contained  in  the  existing  turbulence  is  small.  Therefore  an 
additional  equation  to  describe  the  time  variation  of  v  is  not  needed  nor  do  we  find  a  term  in  the 
energy  equation  to  show  the  heat  generated  by  the  decay  of  turbulence.  No  terms  proportional  to  the 
dynamic  viscosity  appear  because  the  expansion  is  symmetric.  The  heat  conduction,  k,  is  neglected 
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during  the  turbulent  phase  in  comparison  with  the  anomalous  conduction,  pc\a.  We  may  rewrite  these 
equations  in  forms  appropriate  to  a  cylindrically  symmetric  expansion  (hereafter  dropping  the  averaging 
bars). 


*1L 

dr 


+  —  4~  <rpv)  - 
r  dr 


11^ 

r  dr  9r 


0 


(11) 


pc. 


U +  V9I 

dr  dr 


.  DT  1  3  1  \  I  9 

pT  -  Po  T„. 


rpc\  a 


dr 

dr 


(12) 

(13) 


The  approximate  scaling  of  the  channel  expansion  with  time  may  be  inferred  from  these  equations 
by  assuming  the  radial  profiles  of  T—T0,  p—p„ ,  and  v  each  expand  in  a  pseudo-seif-similar  fashion. 
(These  profiles  share  a  common  characteristic  radius,  S,  at  any  time  r  under  the  presumption  of  con¬ 
stant  pressure.)  This  assumption  allows  us  to  replace  the  derivatives  by  simple  proportionalities  when 
we  apply  these  equations  to  a  point  coinciding  with  the  "edge"  of  the  channel.  (See  the  Appendix)  The 
various  constants  of  proportionality  ( a  through  e)  are  retained  when  we  make  these  substitutions  in 
order  that  additive  cancellations  among  terms  may  be  seen  to  be  exact.  The  equations  now  reduce  to  a 
set  of  scaling  relationships 


(p~Po )  ,  Lvp  <P-P«) 

a -  4-  b~f-  -  a  a - 5 - -  0 

t  8  8 2 


pc, 


( T—  T„)  .  ,  <r-  T0) 
c - - - +  d\ - r - 


OT  v  it-  r„) 

+  ep„  R  T„  —  -  cpcv  a  —2 - -  0 


pT-  p0T0. 


(14) 

(15) 

(16) 


The  three  scaling  relationships  would  appear  to  be  insufficient  to  solve  for  the  unknowns  T,  p,  v,  and  8. 
However  after  eliminating  v  and  p  as  we  solve  for  8(r),  ( T—  T0)  divides  out  so  that  the  scaling  is 
independent  of  temperature  as  expected  for  a  self  similar  expansion.  The  solutions  that  may  be  found, 

8!-at  (17) 

v-0,  (18) 

where  8  and  v  are  the  mean  radius  and  radial  velocity,  show  that  the  channel  "diffused  away  at  a  rate 

determined  by  the  anomalous  coefficient.  This  behavior  persists  until  the  turbulence  decays.  Such 
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behavior  may  be  seen  in  Fig.  1  where  the  variation  of  the  square  of  the  radius  with  time  for  an  ohmi- 
cally  heated  channel  in  air  is  shown. 

III.  GROWTH  DUE  TO  HEAT  CONDUCTION 


We  now  consider  the  expansion  of  the  channel  due  to  heat  conduction  after  the  turbulence  has 
decayed.  Dropping  the  anomalous  diffusion  terms  and  retaining  the  thermal  conductivity  the  scaling 
relationships  become 


(T-TJ  ,  ,  (T-T0) 
c - +  dv - - - 


For  the  moment  we  treat  k  as  constant.  Then  a  solution  may  be  found  for  8(r) 


(19) 

v  (T-T„)  n 

+  ep„RT0~  -  ck - -  =  0 

(20) 

PoT0. 

(21) 

pJc^fR) 


P„cp 


(22) 


Thus  the  combined  heat  flow  plus  fluid  response  scales  in  the  same  manner  as  ordinary  heat  flow. 


We  must  now  account  for  the  variation  of  the  heat  conductivity,  k ,  with  temperature.  For  a  gas  k 
is  independent  of  density  but  proportional8  to  vT .  The  channel  behavior  now  depends  on  initial  condi¬ 
tions  as  evidenced  by  the  failure  of  T  to  cancel  out  of  the  scaling  relationships  when  we  assume 

k  -  k-JT.  (23) 

An  approximate  integral  of  the  energy  equation  exists  which  is  independent  of  the  previous 
results  and  which  relates  T  to  8  through  the  initial  conditions.  We  take  a  global  view  and  treat  T and  p 
as  representative  values  for  the  entire  channel.  The  initial  heat  input,  Q ,  must  always  equal  the  work 
done  by  expansion.  W,  plus  the  extra  internal  energy  E.  The  work  done  is  the  pressure  times  the 
volume,  at  density  p„,  needed  to  accomodate  the  gas  expelled  to  beyond  radius  5. 
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W  -  p„ 


(tr  82p„  —  rrS2p) 


P„ 


■■  f)„RT„  7r8 ‘ 


1  - 


The  extra  internal  energy  is 


£  *  pc\ir82  (T—  T„ )  =  p„T„t\iT§‘ 


1  - 


Equating  the  sum  of  £  and  W  to  Q  we  find 


c  2 

°min 

82 


where 


(24) 


(25) 


(26) 


S2  _ k _ 

min  T  ' 

wPo  ^0^0 


(27) 


(As  T  —  oo,  s  —  8mm  which  is  the  radius  at  which  a  channel  of  zero  density  and  infinite  temperature 
would  contain  Q  J/m  and  still  be  at  pressure  P„.)  The  previous  scaling  relations,  including  the  behavior 
of  k  with  T  imply 

*2  - 


Po^p 


Eliminating  £  between  Eqs.  (26)  and  (28)  we  find  the  scaling  is  now 


(28) 


W»T-8-m>— ( 

Po<-p 

where  k„  -  KyfT0. 


Under  normal  circumstances  8  >  Smin  and  the  scaling  reduces  to 

82-  — f 
P  o^p 


(29) 


(30) 


This  limiting  slope  would  appear  as  horizontal  on  Fig.  1*. 


*A  geometric  factor  which  depends  on  the  particular  self-similar  shape  and  choice  of  an  'edge'  has  to  this  point  been  assimilated 
into  a  For  a  Gaussian  profile  with  the  l/e  point  as  the  edge,  equations  (17)  and  <J0>  would  be  of  the  forma*  -  4a r. 
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IV.  APPLICATION  TO  OHMICALLY  HEATED  CHANNELS  IN  THE  ATMOSPHERE 

Figure  1  displays  data  (8,/)  for  an  ohmically  heated  channel  that  has  been  used  at  NRL  to  study 
the  propagation  of  relativistic  electron  beams  in  pre-formed  channels9.  These  channels  are  formed  by 
guiding  an  electric  discharge  from  a  Marx  generator  (U~  250kv,  I  —  \QkA)  with  a  laser-induced,  air 
breakdown  and  have  been  produced  with  lengths  up  to  ~2m.  The  ohmic  deposition  is  ~3  J/cm  and 
the  channels  expand  to  a  ~l  cm  radius  in  ~30ps  before  they  reach  pressure  equilibrium  and  stabilize 
in  size.  After  —  1 50 ps  turbulence  becomes  evident  in  Schlieren  photographs  of  the  channels  and  they 
simultaneously  begin  the  growth  evidenced  in  Fig.  1.  The  kinetic  energy  in  the  turbulence,  ~1CT3 
J/cm,  is  a  small  fraction  of  the  internal  energy  of  the  channel.  It  nevertheless  strongly  influences  the 
channel  behavior.  Rapid  channel  growth  continues  for  —2 ms  with  an  effective  diffusivity  of 
—  500  cm2/ sec  which  is  —  2000  times  the  thermal  diffusivity.* 

The  temperature  of  the  channel  strongly  affects  its  electrical  conductivity  and  thereby  the  propaga¬ 
tion  of  an  electron  beam.  In  these  experiments  the  electron  beam  was  expelled  from  the  channel  at 
!  —  100ps  but  showed  slightly  enhanced  propagation  in  a  channel  at  500 ps  which  demonstrates  the 
importance  of  this  anomalous  cooling  mechanism. 

V.  CONCLUSIONS 

The  treatment  of  turbulence  as  a  diffusive  phenomena  is  often  used  in  engineering  analyses  to 
yield  useful  semi-empirical  results.  Reasonable  experimental  evidence  exists  for  applying  a  similar 
approach  to  the  cooling  of  a  turbulent,  reduced  density  channel.  The  gross  behavior  of  the  channel  is 
thus  "explained" .  It  is  implied  for  example  that  the  radial  density  profile  of  a  Gaussian  channel  with  ini¬ 
tial  central  density,  p,,  and  initial  characteristic  radius,  8,.  evolves  in  time  according  to 

(s  I2  M 

(p-p0)-(p,-po)jy|  e1*1  (31) 

_  7 

*ln  comparing  the  turbulent  expansion  to  true  heat  flow  a  Gaussian  channel  proAle  is  assumed  and  allowance  is  made  for  the  ac¬ 
companying  factor  of  4  in  Equation  <  l7). 
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where 

82  =  5;  +  4a/  (32) 

The  value  of  the  eddy  diffusivity,  a,  remains  however  an  empirical  number.  The  treatment  presented 

here  assumes  a  constant  specific  heat.  c\  .  Such  an  assumption  is  valid  for  channels  in  air  provided 

T  <  3000  °K. 
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VII. APPENDIX:  THE  RADIAL  AMD  TEMPORAL  GRADIENTS  OF  SELF  SIMILAR 
SOLUTIONS 

The  radial  profile  of  any  variable,  K  (eg.  T—  T0,  p -p„,  v,  pv)  is  here  assumed  to  expand  self- 
si  mi  larly 

where  simple  scaling  with  time  is  assumed  for  both  the  overall  amplitude  /  and  the  scale  length  5  such 
that 


/(/)  —  fit",  (34) 

SI/)  -  y/m.  (35) 

and  the  shape  is  always  that  of  g<x).  The  "edge".  x?>  is  a  fixed  value  defined  by 


where  «  is  a  fixed  fraction. 


g<Xe)  *  «  g(0) 


(36) 


Consider  first  the  radial  derivative  of  V,  at  a  point  r,  at  the  moment,  /,  when  the  "edge"  passes  by 

dF 
dr 


8 


}  dx  dr 


hr"''- 


f(')p- 


*<■ 


l 

5  ( / ) 


(37) 
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However  the  slope  of  x(x>  at  the  edge  is  proportional  to  gix,->/Xv 


dg_ 

dx 


fc,  Xt’ 


(38) 


Thus 


d£ 

dr 


j_  fit)  gig?)  =  j_  Hr.t) 


8(/)  x,.  dir) 

i.e.,  the  spatial  derivative  is  proportional  to  the  edge  value  of  V  divided  by  the  radius  ot  the  edge. 


(39) 


Consider  now  the  temporal  derivative  of  V 


d£ 

dr 


f-x..  J  dr 


+  g 


M 

dr 


7  ~  '< 


Starting  with  the  first  term 


J 


dr 


/- 


■  dg  dx 


I*  *  «,  dx  dr 


'  *,• 


fUK^i  & 
X,-  dr 


where  we  have  used  Eq.  (38). 

However  from  Eq.  (35) 

x-T? 

Thus 

dx 

—mr 

dr 

f-x.  =  y'm+l 

7  -  x,  t 

~rwx,- 

l 


The  first  term  thus  becomes 


fdg 
J  dr 


.  Vlr.t) 

-  -/»£ - 

It  -  t 


(40) 


(41) 


(42) 


(43) 


(44) 


We  now  consider  the  second  term.  In  light  of  Eq.  (34) 


& 

*  dr 


7  '  x. 


gpnt 


n-  I 


,  ,  fit)  Vir.t) 

ng{xJ  ”  n - - — 


(45) 


Combining  the  results  for  the  first  and  second  terms  we  find 
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d± 

dr 


I  ■ 


,  V(r.t)  ,  i ■(/-./) 
-ml, - —  +  n - 


Tir.t) 

r 


(46) 


so  that  the  temporal  derivative  is  proportional  to  the  edge  value  of  I-  divided  by  the  elapsed  time  r. 


The  following  proportionalities  are  used  in  the  main  text 


dp  dip  -  p0)  <p-p„) 

dt  dt  J  / 

(47) 

—  ~(rpv)  -  b 
r  dr  6 

(48) 

dr  (T-  Tn) 

dt  ~  1  t 

(49) 

dT  ,  IT  -  To) 
dr~J  8 

(50) 

13/.  v 

—  —  (rv)  =  e  — . 
r  dr  8 

(51) 

The  coefficients  for  the  second  derivatives  are  determined  by  the  requirement  that  the  equations 

reduce  to  the  proper  mass  and  thermal  diffusion  results  in  the  absence  of  the  terms  reflecting  fluid 
motion 


a 


Ii_ 

r  dr 


=  aa 


(P  ~  Po) 
S2 


IJ_ 

r  dr 


dT 

rpc+!7 


cpc\a 


IT  -  T0) 
S2 


The  various  coefficients  are  interrelated  through  the  equation  of  state, 
equation  we  find 


(52) 

(53) 

Taking  derivatives  of  this 


dp  _  ~PqTq  d  T 
dr  ~  T 2  dt 

dp  _  ~PqTq  d  T 
dr  ”  T2  dr' 

Dividing  Eq.  (54)  by  Eq.  (55) 

d£_  d T 

dt  dt 

dp  d  T 

dr  dr 

but  in  terms  of  our  proportionalities 


(54) 

(55) 


(56) 
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so  we  find 


a 


b 


<P  ~  PqJ 
t 

(p  ~  pn) 
8 


c 


( T  -  T») 

t 


,  IT  -  T0) 
d - ; - 


ad  6c 

which  we  employ  in  reducing  the  scaling  relationships. 
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